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ON THE FORM AND STRENGTH OF TREES: 


PART III. THE SECONDARY BRANCHES AND THEIR 
RELATION TO THE PRIMARY BRANCHES 


I. OPATOWSKI 
ARMOUR RESEARCH FOUNDATION, ILLINOIS INSTITUTE OF TECHNOLOGY, CHICAGO 


The secondary branches are considered in this note mainly for the 
ln of applying the principle of maximum strength to the primary 
ranches 


1. A reason for secondary branches. If the secondary branches 
are treated as horizontal cantilevers, they may be subjected to the 
same mathematical analysis that was carried out previously (Opa- 
towski, 1944b) for the primary branches. The same function K rep- 
resents here the advantage in mechanical strength of a primary 
branch with many secondary branches with respect to a biologically 
equivalent primary branch with only one secondary branch. 

2. A branch of maximum strength. The principle of maximum 
strength, which was applied to the trunk, (Opatowski, 1944a) may 
be used in a similar manner for the primary branches. It must be 
supplemented by the condition of maximum length of the branch, in 
the same way as in the case of the trunk the condition of maximum 
height had to be added to obtain the shape of the trunk. The condi- 
tions of maximum height for the trunk, and of maximum length for 
the branch, may both be justified by biological considerations. 

We consider the primary branch as a horizontal cantilever which 
is bent under the action of its own weight and of the load due to sec- 


Fos 


strength of materials. Since we do not discuss the growth of the 
branch, we may consider its weight as a fixed and known quantity. 
Then the conditions of maximum strength and of maximum length 
require that the maximum bending stress of each cross-section be the 
same throughout the cantilever, that is, that the latter be of uniform 
strength (see for instance, Timoshenko, 1940, pp. 209-210). 

One may think to take into account the weight of the secondary 
branches by a suitable increase of the specific weight of the primary 
branch. However, the ordinary theory of strength of materials would 
be in such case inadequate to deal with the problem of uniform 


1 


case (12308 3/19/42 Ad v7 


UC 


ondary branches. We apply the ordinary methods of the theory of - 
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strength, because it would lead to the inacceptable conclusion that 
the deflection at the free end of the branch is either imaginary or is 
infinity (Opatowski, 1945). Therefore, the secondary branches must 
be considered as concentrated weights. A similar situation arose in 
discussing a trunk of maximum strength (Opatowski, 1944a). 

Let us number the secondary branches progressively: 1, 2,--- , 
i,---, n, starting from the free end. Assume that the cross-section 
of the primary branch is circular and that the action of each second- 
ary branch is equivalent to a concentrated weight Bb"; acting at a cer- 
tain point P;. Let B’; be the weight of that part of the primary 
branch, alone, which joins the points P;., with P;, and J; its length, 
that is, the distance (Pi,,, Pi). All the weights B’; and B"; are as- 
sumed to act in one vertical plane. Call B; the resultant of all weights 
acting between and at P,, P,,---, P; and M; the resultant moment of 
all those weights with respect to P;. 

Instead of weights and moments we will use, for simplicity, the 
following quantities which have the dimensions of areas and volumes 
respectively: 


Bi =FB',, Bb’, =FB";, Br= FB, 4—=FM,, 
where [Opatowski, 1944a, equation (3) ], 
F=4f'/(ao,) = 103 X 10° f'G*. (1) 


We will maintain the name of forces for /’s and of moments for 
ws. The units are throughout the paper kg and em. We have the fol- 
lowing obvious relation: 


Bir = Bi + Bit Bin. (2) 
The weight acting at the free end of the primary branch is B”, and 
we must assume that it is # 0 , otherwise we would run again into an 
infinite or imaginary deflection. This detail however may be easily 
taken care of by considering the very tip element of the primary 
branch as a secondary branch. 
If r; is the radius of the cross-section of the primary branch at 
P; , the condition that the maximum bending stress in that cross-sec- 
tion be equal to «,/f’ is (Timoshenko, 1940, pp. 209-210) 


w=, (3) 


whereas at a point P between P; and P;., , at a distance x from P;, 
that condition is (Blasius, 1914; Opatowski, 1945) 


mt Paty] (@- OO @LAs=[rayk, (4) 
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where r(x) is the radius of the cross-section at # , and 
y =m f'G/ (250 o) = 8238 X 10° m'f’'G- . 


Equation (4) is an integral equation in r(x), which differentiated 
twice with respect to x gives 


srr’ + 6r2 yr, (5) 


where the primes stand for derivatives with respect to x. The initial 
conditions are: 


r(0) = 15, 7 7°(0) S6i/ (878): (6) 


The solution of equation (5) under the conditions (6) is the hyper- 
elliptic integral 


2=3 [ Zoonpr dr, (7) 


where 
eS Ty) = [1.24 (2° — 75) + Biz]? . 


' The weight of the part of the primary branch between P; and P,., is 
easily obtained from equation (7): 
Bs=Z (Tin, Ti) — Bi. (8) 

The equations (2), (8), (7), and (8) relate the radius of the 
cross-section of the primary branch to the weights B"; of the second- 
ary branches and to their mutual distances 1; . These equations provide 
a possibility of testing the principle of maximum strength on an actual 
branch and may be used also for calculating the safety factor f’. The 
applicability of the formulae is not limited to horizontal primary 
branches, because they may be used also for horizontal parts of curved 
branches. The theory could be easily extended to inclined or curved 
branches, but the purely statical considerations of mechanical strength 
are not sufficient to explain the curvature of a branch. 

The function r(z) is completely determined by the quantities 
B",, l; and all the 7;’s may be calculated successively from the equa- 
tions (2), (3), (7), and (8). Such calculations involve however an 
evaluation of the hyperelliptic integral (7) which may be carried out 
through a power series expansion of r(x). We obtain directly from 
equations (5) and (6): 

Pe) =Tet Fe 9/2) 22 + (7% 1/6) 3 +. (9) 
where 
Y; = p;/ (37,7) ; ee, = (y/3) a (27;2/7;) ; 
r= (y — 169’ 1) 7. / Br): 
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According to the theory, the secondary branches do not cause discon- 
tinuity in the radius of the cross-section, because the bending mo- 
ment increases continuously from the free end towards the junction 
with the trunk (equation 3). One gets vi from equation (9) putting 
“Lx L; - 

As a numerical example, let us consider a horizontal piece (p) 
of a primary branch between two consecutive secondary branches. 
Cut off the primary branch at a point P; of (p) and let 7; = 5 cm be 
the radius of the cross-section at P;. Let B; = 5 kg be the weight of 
the part of the branch cut off (complete, that is, with secondary 
branches and leaves) and M; = 600 kg cm the moment of B; with re- 
spect to P;. Let m’ = 2.5, G = 0.5 so that G/* = 0.84. From equa- 
tions (1) and (8) we get the safety factor 


f' = 972 G4 7r2/M,; = 85. 


Therefore, y = 0.0008181 , F = 0.2083, 6; = 1.042. For 1; = 25 cm 
we get from equation (9) 


’ 


Tie = 5 + 0.347 + 0.061 — 0.005 + ---=5.4 cm. 


The convergence is good. 
Since y is small one could expand expression (7) in power series 
of y, which would give 


x = 6; (e — r;°) + pr (0.67;° rs — (2.22573 — 0.3757;5) y 5 ae 


However, the convergence of this series is not better than that of ex- 
pression (9) and the high powers of 7 make its use inconvenient. 

The writer is indebted to Professor N. Rashevsky for his contin- 
uous interest in this work, to Professor A. S. Householder for the op- 
portunity to discuss the main ideas of this paper at a seminar meet- 
ing of the University of Chicago, to Dr. H. D. Landahl, Mr. J. Nielsen, 
and Dr. W. 8. McCulloch for valuable criticism and suggestions. 
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MUSCULAR DYNAMICS AND MUSCULAR EFFICIENCY: 


I. THE ISOMETRIC LENGTH-TENSION DIAGRAM 
OF STRIATED SKELETAL MUSCLE 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


The isometric length-tension diagram for individual fibers and for 
whole muscle is considered, and it is proposed that the tension » may be 
represented for any muscle whose ‘fibers are parallel and not in series, 
in the form 

p—f(x)+ B¢(a,l,x), 


where the form of ¢ is known and the same for all muscles, the para- 
meters a and / are experimentally measurable, ~« represents the length, 
B the degree of activity of the muscle, and the form of f(x) varies from 
muscle to muscle. 


We wish to consider the muscle as the source of power for the 
operation of the bodily machinery. Consequently, our chief interest 
shall lie in the muscle as a whole, in situ. The importance of muscle 
lies not only in the fact that it alone drives the machinery in the per- 
formance of overt activity. On the contrary, even if we do not go so 
far as M. F. Washburn (1916) and identify all mental imagery with 
muscular movement, the work of E. Jacobson (1934) and others has 
made it abundantly clear that muscular activity of at least an incipi- 
ent sort is intimately associated with even the “higher” thought pro- 
cesses, while at the same time certain fairly successful theories of 
discriminatory processes (Rashevsky, 1940; Householder, 1940; 
Householder and Landahl, 1944) have postulated for muscles a role 
very similar to that which M. F. Washburn attributes to it. 

In nerve, considered with respect to its role in the reactions of the 
organism, the interest centers about the action-spike, its conduction 
and its transmission. Because of the all-or-none law and the phe- 
nomenon of refractoriness, nervous activity is quantized (McCulloch 
and Pitts, 1943), the principles of nervous coordination are express- 
ible in terms of merely the occurrence or non-occurrence of each of a 
certain class of events (the firing of specified neurons) in each of a 
class of time-intervals (covering the life-time of the organism). It is 
true that the complexity of the structure of the nervous system and 
the shortness of the elementary time-intervals must necessitate the 
employment of statistical analyses (Landahl, McCulloch, and Pitts, 
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1943; Householder and Landahl, 1944), but the bases for a theory of 
nervous coordination seem to be conceptually simple. 

But for the present, at least, there exists no such reduction of 
muscular activity to simple elements, and the great divergences of 
opinion among investigators, and their sharply discordant results, 
make it very difficult to discern a direction in which to look for a 
reduction. 

We wish to consider the mechanical response of muscle, having 
in view principally the whole muscle in situ. The main object will be 
to arrive at quantitative formulations for describing, first, this re- 
sponse as a function of time and a functional of the stimulus, and 
second, the efficiency of the muscle, regarded as a machine for con- 
verting chemical energy into mechanical work, under the various con- 
ditions to which it may be subject. With this purpose in view we 
shall first examine, more or less in isolation one from the others, sev- 
eral sets of properties of the muscle, describing such formulations as 
have been proposed in the past, and attempting to revise or supple- 
ment these where possible. 

The simplest set of properties is contained in the classical iso- 
metric length-tension diagram (cf., e.g., Fenn, 1925; Hill, 1926; Ram- 
sey, 1944). In the diagram for resting muscle, the muscle is stretched 
by a given amount, and the resulting tension is plotted as ordinate 
against the length as abscissa. For the active muscle the muscle is 
first stretched, then subjected to tetanic stimulation without being 
permitted to shorten, and the maximal tension is then recorded. Upon 
stimulation the tension is ordinarily raised, but according to E. Weber 
(“Weber’s paradox’, 1846; reported, e.g., by W. O. Fenn, 1925) when 
the stretch is sufficiently great the tension may diminish. 

Whether it is a single fiber or the whole muscle under considera- 
tion there must result a set of curves which depend upon at least one 
parameter for characterizing the stimulation (cf. Th. Christen, 1911). 
Let us consider one of them, supposing the temperature and other 
conditions to be held constant, and writing its equation 


p—f(x). (1) 


There are other ways in which one could obtain a length-tension 
diagram. One could first stimulate and then stretch to length x , and 
one could vary the rate at which the stretch takes place, etc., and 
there is no reason a priori for supposing that one would always come 
out with the same tension p at any given length x, even though the 
muscle may at no time have been stretched enough to injure it. Thus 
EK. Asmussen (1936) finds that if a fiber bundle is first stimulated and 
then stretched beyond about 25% of its normal length, without stimu- 
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lation and without tension, the tension becomes increasingly greater 
than that which is developed when the stimulation follows the stretch. 
This effect he associates with the fact that he finds the threshold of 
stretched muscle to increase with the amount of stretch, so that we 
may infer that in the prestretched bundle there are fewer fibers acti- 
vated by the stimulus, and hence less tension developed. But it is not 
clear from E. Asmussen’s results that, once a fiber is activated the 
subsequent dependence of its tension upon its length is in any way 
dependent upon the initial length or that hysteresis is otherwise evi- 
dent to any substantial degree, nor is this evident from the work of 
previous investigators when the stretches are not extreme (Fenn, 
1925). Moreover, while R. W. Ramsey and 8. F. Street (1941) find 
that when a fiber is allowed to shorten too much (below about 2/3 its 
length at zero tension without stimulation) it then enters irreversibly 
a new condition (the “d-state’’), they deny the presence of hysteresis 
where the changes are not extreme. D. Brown (1936) reports, indeed, 
that the tension developed in isometric twitch increases with the 
initial tension, but it does not follow that this is the case in tetanic 
stimulation. It seems reasonable to suppose, therefore, that if we 
exclude extreme degrees of stretch and shortening (and perhaps ex- 
treme rates of stretch and shortening) the hysteresis will be negli- 
gible, and we may regard the unstimulated or the tetanically stimu- 
lated muscle or fiber as an elastic body for which the tension is a 
single-valued function of the length as given by equation (1). 

If our supposition is legitimate, we may define the potential 
energy of the muscle at length x , when subject to the given mode of 
sustained stimulation, as 


w (2) = [ teas, 


where x, is the lower limit of the range of reversibility. Alternative- 
ly, we may extend the range of definition of f(x) for x < a and de- 
fine f(x) on this range as the tension developed when the muscle is 
stimulated at the length of zero tension and then allowed to shorten 
to the length «. Then we may replace the above function by 


w(a)= f seas. (2) 


Clearly, then, if the muscle is allowed to shorten from a length 
x, toa length x, it is capable of doing work in amount w(x) — w Cas 
If the muscle is allowed to shorten isotonically under tension f (x2) 
it will do mechanical work amounting to only (#1 — 22) f (#2), whence 
the difference must be dissipated as heat. Since w’ = f/f, this differ- 
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ence is, by a standard theorem, 
w(t) — W(t) — (ar — He) f (%2) = F(X — Xe)? F(X), (3) 


for some « satisfying 2, = « = %. The amount of heat per unit of 
shortening is therefore 4(%, — %2)f'(x), a quantity which vanishes 
with the amount of shortening. If the shortening takes place under 
some lesser, fixed, tension p and is suddenly arrested at the length 
2, then the heat liberated per unit of shortening is 


f (te) —p + 4(%1 — Xe) f' (22). 


This we may denote the viscous heat, though it is quite independent 
of the measure of the viscosity. 

Now A. V. Hill (1938) finds that the amount of heat which is 
actually liberated by the stimulated muscle in the process of shorten- 
ing—over and above the “maintenance heat” which is developed while 
no shortening occurs—is independent of the work which is performed, 
which is certainly not true of the viscous heat just defined. Granting 
this result, we must seek the explanation in the molecular mechanism 
of contraction (Brown, 1941), and in the ordinary thermodynamic 
effects (Hill and Hartree, 1920). 

Determinations of the length-tension diagram have been under- 
taken by a number of investigators (see W. O. Fenn, 1925, for the 
earlier papers). Measurements by E. Weber (1846; reported by W. 
O. Fenn, 1925; Th. Christen, 1911; and others) indicated that 


d*p/dx? > 0: (4) 


the increase in tension per unit of stretch increased with the amount 
of stretch (“Weber’s law’). Th. Christen (1911) presents alterna- 
tive theoretical functions for whole muscle, both satisfying (4), and 
both based upon formulas proposed by previous authors, but modi- 
fied to eliminate certain formal difficulties. The first of these is taken 
from G. Weiss. Th. Christen’s form is 


du/dp=a/[p+e], p=c(e*/4—]), (5) 
This is equivalent to 
dp/dx = B e/*, (6) 


according to which the ratio of the increment of tension to the in- 
crement of length is an exponential function of the length. 
Th. Christen’s other equation has the form 


dp/dx=c + a(x—s), p=c(x%—s) + fa(x —s)?2, Cy 


where the ratio of the increments is linear instead of exponential. 
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M. G. Banus and A. M. Zetlin (1938), on the basis of extensive 
measurements upon unexcited whole muscle (cat and frog gastrocne- 
mius, frog sartorius, cat rectus), propose, in present notation, 


log p=kz+a, 


and find very close agreement. But this is essentially G. Weiss’s equa- 
tion which Th. Christen replaces by equation (5) in order to eliminate 
the absurdity at p = 0, and if we make a proper choice of the para- 
meters, we may do likewise. Now M. G. Banus and A. M. Zetlin find 
that if the muscle-fibers are removed from the connective tissue 
sheath, a curve of almost identically the same form is obtained from 
this sheath alone. This may indicate either that virtually the whole 
tension of the resting muscle is due to the connective tissue sheath 
alone, or else that the length-tension diagram of the individual un- 
stimulated muscle-fibers has the same form as that of the connective 
tissue sheath. While the somewhat earlier work of F. Sichel (1934) 
seemed to indicate that the diagram for the individual unstimulated 
fibers is linear, E. Asmussen (1936) and R. W. Ramsey and S. F. 
Street (1940) criticize his methods and obtain, for individual fibers 
and for small bundles, curves which closely resemble those of M. G. 
Banus and A. M. Zetlin. According to R. W. Ramsey and S. F. Street, 
while some of the curves they obtained were exponential they were 
not all so, although they did all exhibit the upward concavity in ac- 
cordance with Weber’s law. 

R. W. Ramsey and S. F. Street (1941) attribute all the tension 
exerted by the unstimulated stretched fiber to the sarcolemma, since 
they have succeeded in severing the myo-fibrils by local pressure, thus 
obtaining a segment of clear sarcoplasm, and found the length-tension 
diagram of the injured fiber to be almost identical with that obtained 
for it prior to the injury. A further detail which strongly supports 
this supposition is the fact that parasites have been observed freely 
moving within the fiber (Cowdry, 1938). Both observations seem to 
throw doubt on the measurements of F. Buchthal and J. Lindhard 
(1938-39) indicating that during passive stretch the isotropic and 
the anisotropic segments are stretched in different proportions, unless 
we are to suppose that the sarcolemma has itself a segmental char- 
acter. 

A. Bairati (1937) finds by histological examination that the sar- 
colemma consists of virtually inextensible anastomosing fibrils im- 
bedded in a matrix, these fibrils running, some parallel to the axis of 
the fiber, some obliquely. His conception seems to be that the fibrils 
which run obliquely form a kind of network with meshes which we 
may picture schematically as rhombic, that the longitudinal fibrils are 
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attached to opposite vertices of these rhombi, some being connected 
to one, some to the other tendinous termination of the fiber, and that 
the tension in stretch is due to the distortion of this net. Before 
stretch these fibrils are described as wavy. A. Bairati and also R. W. 
Ramsey and S. F. Street (1940) suggest that the form of the length- 
tension diagram of the resting fiber may be explainable on the basis 
of the gradual taking up of the slack and subsequent stretch of these 
fibrils. However, some of these fibrils must be taut even at the rest- 
ing length since otherwise we could not understand the return to this 
length (on this basis) when the external tension is released. We can- 
not introduce sliding friction, as A. Bairati seems to do, to explain 
the tension in stretch for an analogous reason—some elastic tension 
is required to return the fiber to its resting length. Now if we sup- 
pose that the longitudinal fibers are inextensible and only the oblique 
fibers are extended, since the length of an oblique fibril is of the or- 
der of the radius of the fiber, this being only about one-tenth the 
length of the fiber and of the order of one-tenth of the amount of 
stretch the fiber can sustain reversibly, the extent to which these 
oblique fibrils would be stretched would be quite enormous. It seems 
much more likely that the fibrils act as a safety catch designed to 
prevent altogether a stretch exceeding.the elastic limits of some other 
elastic mechanism which is normally operating. 

This leaves us, by successive elimination, with the matrix in 
which the sarcolemmal fibrils are imbedded as the only possible elas- 
tic mechanism in the resting muscle-fiber. Apparently but little is 
known of this, and there seems to be little profit in presenting any- 
thing but a phenomenological equation such as Th. Christen’s for de- 
scribing its elastic behavior. But there is certainly no occasion for 
surprise in the failure of Hooke’s law over the wide ranges of stretch 
employed, and an explanation of its failure is neither more nor less 
essential then would be an explanation of its validity if it had been 
found to hold. It is, however, interesting to note the resemblance of 
the length-tension diagram to that of rubber (James and Guth, 1944), 
especially in view of the well-known fact that both muscle and rubber 
contract on heating (Hill and Hartree, 1920). 

When the muscle-fiber is stimulated the length-tension diagram 
differs quite considerably, passing through a maximum and a mini- 
mum and later following the resting diagram quite closely. While EB. 
Asmussen (1934, 1936) and others think that the all-or-none law is 
a property of the motor end-plate but not of the fiber, the greater 
number of investigators seem to be of the contrary opinion (Gelfan 
and Gerard, 1930; Ramsey and Street, 1940), and we shall adopt the 
latter hypothesis. Hence associated with any fiber are but two dia- 
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grams, that corresponding to the state of rest and that corresponding 
to the state of complete activity, at least if the frequency of stimula- 
tion is high enough to prevent relaxation between stimuli. According 
to R. W. Ramsey and S. F. Street (1940, 1941) when the ordinates of 
the diagram for rest are subtracted from those for activity, and the 
ordinates and abscissas taken to represent fractions of maximal ten- 
sion and of rest length, respectively, the resulting differences agree 
very closely from fiber to fiber, all possessing a maximum at (1, 1) 
and intercepting the horizontai axis at 2. According to R. W. Ramsey 
(1944) the failure of this curve to pass through the origin is due 
chiefly to the resistance of the sarcolemma. 

Presumably, therefore, the isometric length-tension diagram of 
the whole muscle should involve a single parameter only, which can 
be taken as the proportion of fibers simultaneously activated, contrary 
to Th. Christen’s (1911) assumption that two are required. It is hard 
to see how two independent parameters could be involved unless we 
suppose the distinct fibers to have sensibly different curves and to 
be provided with a neural mechanism capable of stimulating the fibers 
independently. Th. Christen argues that one can voluntarily exert a 
predetermined (limited) amount of tension irrespective of position 
and one can also assume a predetermined position (with muscles at 
fixed lengths) irrespective of tension. But for this to be possible one 
need only suppose that two types of kinaesthetic receptors are avail- 
able, one for tension and one for position, both suitably linked with 
the motor neurons—the single parameter may itself be a function of 
two others—and there is evidence for such distinct receptors (Mat- 
thews, 1933). 

If the contractile mechanism is linked in parallel with the static 
elements which operate alone when the muscle is at rest, if the muscle 
consists of parallel fibers (none in series) of approximately equal 
lengths, and if the tendinous connections are of very high elasticity, 
let us write for the whole muscle at rest 


p=f(x), (8) 
where f is perhaps exponential in form. Suppose that each fiber is 


only a times as long as the whole muscle, the rest consisting of ten- 
dinous connections, and let 


a= (&) (9) 
represent the tension taken up by the contractile elements of an indi- 
vidual active fiber, where a represents a fraction of the maximum for 


the fiber and é represents the multiple of the resting length. If / is 
the resting length of the whole muscle, then a / is the resting length 
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of the fiber and therefore the length x of the whole muscle is given by 


Edit ly Mi ot (2 cto 8 
or 
é=[x#— (1—a)l]/(al). (10) 
Hence 


n= ¢{[z— (1—a)l]/(a)}. (11) 


Therefore, if 6 represents the intensity of activation, being equal to 
the number of fibers activated multiplied by the maximal tension con- 
tributed by the contractile elements of each, the isometric length-ten- 
sion diagram of the whole muscle should have the form 


p=f(«) +6 el{le— (1 —a)l]/(al}. (12) 


The parameters / and a are directly measurable; according to R. W. 
Ramsey (1944) the function ¢ is very nearly the same for all fibers. 
There remains only the single parameter # and the additive function 
f(a“) which is independently determinable and appears to be approxi- 
mately exponential in many cases. Note that we here employ no as- 
sumption whatever as to the source (whether sarcolemma or sarco- 
plasm) of the elastic properties of the resting muscle described by 
the function f(x). 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Following the general form for the differential equation of organ- 
ism and colonial growth, there is derived a rational formulation for the 
growth of a bounded cell community (e.g., an organ) equipped with a 
food supply and a waste removal mechanism. It is shown how, from the 
integral form and an empirical curve, the vital coefficients of the equa- 
tion can be derived. Changes to be expected in these coefficients are 
discussed, and the analytic methods for assessing them are set forth. It 
is hoped that these equations and similar ones will make it possible to 
relate empirical curves to the mathematico-biophysical theory of the cell. 


In the following analysis we will examine the growth of a cell 
community, spatially bounded and equipped with a mechanism which 
supplies its food and removes its wastes. In the course of our deriva- 
tions it will be necessary to make many gross assumptions and ap- 
proximations. We feel that no apology is necessary for this policy 
provided that the temporary nature of the work is fully realized. Our 
main thesis will be in line with the notion set forth in a previous pa- 
per (Morales and Shock, 1942): A cell has its own fundamental 
growth rate. In the proximity of other cells it suffers interactions 
which change this rate. The reduction of food supply per cell, irre- 
spective of the mechanism, and the increased concentration of waste 
products, irrespective of the mechanism, always give rise to reduc- 
tions in the growth rate of the community. If one identifies “food” 
with “reactant”, and “waste” with “product”, then the problem is 
reduced to one of the chemical kinetics in which diffusion is recog- 
nized as a factor. In this sense our analysis is a gross version of the 
fundamental theory of cell division developed in Mathematical Bio- 
physics (Rashevsky, 1938; 1940); in fact we hope that analyses of 
the type to be presented will provisionally link experimental work 
with the mechanical theory of the cell, and thus open to the latter a 

* The opinions or assertions contained herein are the private ones of the 


writers, and are not to be construed as official or reflecting the views of the Navy 
Department or the Naval Service at large. 
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vast array of guiding data in the form of growth curves of com- 
munities. 

For the sake of concreteness we will use the terminology proper 
to a well defined homogeneous cell community, such as a vertebrate 
organ, equipped with an arterial and venous food supply. It will be 
clear as we proceed, however, that this interpretation of the situation 
is not unique. Furthermore, we will carry out the derivations con- 
sidering only the interactions of food supply; later we shall readily 
show that the same remarks apply mutatis mutandis to the interac- 
tions of waste accumulation. 

Let N be the number of cells in the organ at time, ¢. Then the 
rate of proliferation per cell is dN/Ndt. Under the assumption that 
there is a fundamental proliferation rate which is subject to the in- 
fluences, for example, of hormones or radiations on the mitotic mech- 
anism, but not to other causes, we have provisionally that 


1 dN 
——- Constant + ---=a+---. (1) 
N dt 
Suppose, furthermore, that the cells ingest from their medium an 
amount of food which is appreciable in comparison with the amount 
supplied to each cell by the blood*. If 


R li. sec. = rate of blood flow into the organ, and 
x” moles li-1=concentration of the i-th nutrient in the arterial 
blood, 
then 
Rx moles sec. = T = total influx of food, and 
T“)/N moles sec.“cell“ = food available per cell. 


Let us imagine that each cell actually ingests B‘® moles sec.cell of 
the 1-th nutrient. Then the total amount of food made unavailable by 
this mechanism on the average is B“ N moles sec.". The resulting drop 
in the steady state concentration of the nutrient must be proportional 
to BN, and therefore the negative effect on the rate of prolifera- 
tion per cell is likewise proportional to it. Since the numerical value 
of B is not to be discussed, we omit the introduction of a proportion- 
ality constant, and write merely 


(Effect)  =—BWN , 
; Next we shall study blockage of nutrient streams by regarding 
it as an immobilization effect which all cells exert on their neighbors 


*In the author’s opinion this could be ascertained experimentally from 


orig: studies such as have been made for certain amphibian organs (Morales, 
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even though they do not actually ingest such food in competition. Thus 
we may consider any nutrient molecule to be in the force field of many 
cells. Consider any two such competitors. The net force with which 
either attracts the molecule is less than it would be were the other 
competitor not present. The net effect is a diminished food supply 
per second for either of the two cells. The same result would obtain 
if, through the influence of the neighbor, the concentration near one 
cell had been lowered by a fictitious amount. Treating the problem 
this way circumvents the more exact but almost hopeless vectorial 
analysis of the situation. 

Choosing spherical coordinates, r, 6, ¢, with origin at the n-th 
cell, we proceed by letting K,‘* (7) moles sec. cell-? be the immobili- 
zation effect of the n-th cell on any cell. Let p(7, 6, ¢) be the cell 
density in cells cm. Then the effect of the n-th cell on all cells of 
the organ will be 


{ K,, pr dv moles sec.~! cell-, 


organ 


and the total interaction effect would be 


n=N 
E> i} Kn pn dd. (2) 


n=1 organ 


However, it is clear that to evaluate each integral the origin would 
have to be shifted from cell to cell for all N cells, and in each case the 
equation of the surface which gives the organ shape (hence the limits 
of integration) will be different. We avoid this difficulty in the fol- 
lowing approximate way: Suppose that for the n-th cell the equation 
of the organ surface is available in the form 


Rn =1n(0,¢)- 


We then define the mean radius of the surface “viewed” from the n-th 
cell as 
25 - £ 
eran J r, 6 dd 
Tn —i ee n ’ 
An f 
0 0 


and the mean of these means as 
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Thus as far as the n-th cell is concerned, we regard the organ as a 
sphere of radius 7) with center on the n-th cell, and the integral (2) 
can be written as ~ 


%, CT 20 


E =3 {| f | K, (7) pa(r, 6, 6) 7? sin 6 drdé d¢ . (3) 


We will regard equation (3) as the general form for the nutrient 


immobilization effect. In the event that p be a constant, say p, inde- 
pendent of coordinates, equation (3) reduces to 


—_ nN As 
E=Anp (i K, (r) 72 dr. (4) 
n=1 


This is certainly a justifiable assumption for many organs. In this 
case, if V is the volume of the organ, the gross density is merely 


p=N/V. 


Although many organs are quite uniform histologically, it is well 
known that at any one time there are important differences in amount 
of metabolic activity from region to region within them, consequently 
K, (r) will generally be different for each cell. If we assume a mean 
value for this quantity, or if for some reason we believe the organ in 
question to be physiologically homogeneous,* then equation (4) may 
be simplified further to 


B= N-{ = ii K® (r) r? dr ) : (5) 


The factor in parentheses includes a definite integral not containing 
the dimension “cell”, consequently it follows that no matter what is 
the quantitative law of interaction, K), the total interactive effect, E , 
ts proportional to the square of the number of celis, and the inter- 
active coefficient or “vital coefficient” (Kostitzin, 1939) is a function 
of the organ or colony size, not necessarily of cell number. The inte- 
gral of equation (5) can be evaluated for some simple laws of action, 
and it is interesting to consider two of these because they may have 
some bearing on the more fundamental theory referred to initially. If 


* This is very likely true of gland i imi 

: IS 18 : glandular tissues. To some extent a similar 
qualification is necessary in considering “B”, but for that case the mean value of 
B is easy to compute from arteriovenous differences. 
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K® (7) — a™ rm (m =: 0), 
4naa | Toe™ 
V 


in which it can be easily shown that m < 3 must hold. If, on the 
other hand, 


po=x| | 2 


om 


K® (7) ——(.) err 
with 
T>0F 
then 


ys = Wise) Awe 5 AM) 
E™ = N? (42a) ——_ ————__—_——__, (7) 
V 
In both equations (6) and (7) it is probable that E will increase 
with time unless the organ is characterized by anomalously large cells 
and a low proliferation rate.* The vital coefficient of N? which we 
may henceforth call C’ such that EH = — C“N?, will tend to di- 
minish in both equations (6) and (7), if we regard V as of the order 
of 7°,. 
In any event, we now have reasonable grounds for setting the 
sum total of this effect as equal to 


2 ioe 275 2 ero 


—C‘) N? moles sec."! cell? . 
Returning now to our original problem, we see that the i-th food 
not unavailable per cell per second is 


Ts) 
— BY N — CN? moles sec. cell“! 


and for f foods ((=—1,2,---f) itis 

i=f TW) 

> — — BON — CON?, 

i=1 

Let us now see to what extent these remarks are applicable to the 

elimination of wastes. Since blood is incompressible, the rate at which 
blood flows out is R li. sec. (actually this should be blood plus lymph, 
but for purposes of formulation this is immaterial). If 7” is the con- 
centration of the j-th waste in the venous blood, then TT” = Ra“) 
moles sec. is the total efflux of waste 7, and T/N moles sec.*? cell 


* This follows because N is of the order of r°,. 
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is the waste removed per second per cell. Suppose each cell produces 
per second an amount of waste equal to B®, The total of such effects 
is BON. For converse reasons to those outlined above, the immobi- 
lization of wastes gives an effect proportional to C’)N?. Now the to- 
tal amount of waste accumulating or affecting the steady state waste 
concentration should be proportional to 


ey OD 
— > Sens + BON + CN? for w wastes ,1,---w. 
jet 


The rate prescribed by the tentative equation (1) is now modified by 
the positive effect of excess food, and the negative effect of excess 
waste to give 


aN balay oa SZ. 
et ee ee — BYON — C#N? 
dae = ( N C 
j=w TO) (8) 
-> le: ee PUA cons) : 
j=l N 
or, rearranging and using the fact that T= Rz, 
Seek ee 
aN é j 
—— a+ Rk ——————— 
dt N 
(9) 


= > Bw + = B®)N —(SCH+43 onse| m 


To our degree of approximation we are not justified in attempt- 
ing to distinguish between the effects of various metabolites, as would 
be suggested by the form (9). Accordingly we substitute each of the 
expressions in parentheses of form (9) by gross coefficients, A, B, 
and C , remembering, however, the general nature of these quantities. 
This leads finally to the differential equation of the growth of the 
organ: 

aN 
Pienaar eS (10) 


Before such an expression can be of any but speculative value, meth- 
ods must be developed to test it experimentally. In the case of equa- 
tion (10) this can fortunately be done in complete fashion. As we 
shall discuss below, variations in A, B, and C are to be expected in 


~ 
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the course of the life-span, but in an approximate general way (or 
under certain experimental conditions, in a rigid way) the general 
integral of equation (10) will not differ seriously from the integral of 
(10) in the special case of A , B, and C are constants. As a first step, 


then, we will study equation (10) under the assumption of constant 
vital coefficients. 


The slope at the origin is evidently zero, as it should be for ana- 
lytic reasons (Winsor, 1934). It is again zero when 


—B + B+ 4AC 
N=—____—_ 
2C 


Of these only the positive root has physical meaning. We denote this 
positive root by N.,. Differentiating equation (10) once more, 


d?>N 
dt? 


Evidently d?N/dt? = 0 when dN/dt = 0, i.e., at the origin and for 
N =-N,,. It is also zero when 
—B + VB? + 3AC 
3C 
We denote the only physically significant root by N;. Thus it appears 
that similar to the simpler so-called “logistic” function, this equation 
gives one asymptote and one inflection between zero and infinity which 
need not be so that 2N; = N,,. 


The expressions for N., and N; make it possible to write B and 
C in terms of these quantities and of A. Thus 


1 2N,—-Na 
Ni Nw 2Na — 3N1 
1 Ne ONG 
——_——_A 
NiNw 2Na eS 3N, 


Let us now consider the direct integration of expression (10). This 
yields by elementary methods 

N? B 2CN-— 6 
og —————__ + —— ‘tanh —_—=At+@G@, (12) 
A—BN-CN? vV-@ Nf et | 


where G’ is a constant of integration and —qg = B? + 4AC. Making 


=— (8CN? + 2BN — A). 
dt 


(11) 


— 


4] 
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use of the expressions derived above for Noo and N;, we readily trans- 
form expression (12) to give 
N? 
0g SS. eee Bh a 
NoNr(2Na 38ND) — (Na? — SNP )N Nr Na 


Nw” a 3N/? 
ae - 
(New — 3Ni) (Nx. — Nx) 
2(2N, —s N.)N == (NS ae 3N 77) 
NoMa ee 
(383N; — No) (No — Nz) 


It follows that the constant A and the intercept G are obtainable from 
a plot of the appropriate function of N,, N.., and N vs. time; conse- 
quently, in view of expression (11), all vital coefficients are deter- 
minate from experimental data. This concludes our derivation for 
the special case. 

When, in accordance with expectation, A , B, and C are varying 
in time, their changes can be followed by analyzing the empirical 
curve. A method of doing this for the general growth equation is 
discussed elsewhere (Morales, 1944, in press) ; so it will suffice here 
to present the application to the particular problem at hand. If the 
differential operator, H , is defined as 


tanh" (13) 


=At+G@. 


Pe gh | 
yi Lad 
dN dt 
dt 
then, for instance, at (NV, t): 
aN 
N H° — N? 
dt 
aN 
1 Hi — 3N2 
dt 
dN 
‘tee H? — 6N 
dt 
BS ee SS 
N N? N38 
1 2N 3N?2 
0 2 6N 


and analogously for A and C. Thus it is clear that from the em- 
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pirical data it is always possible to obtain the vital coefficients. 


Let us now return to a consideration of the physiological basis 
of our equation: 


= ( o-57 )-BN— on |. 


If our hypotheses hold exactly, then of course relation (12) will hold 
exactly. However, this is not to be expected, nor in fact is it particu- 
larly desirable in the sense that a perfect fit would not substantiate 
our analysis as well as a fit which deviated in a certain way. The 
reason for this lies in the fact that we have not presumed the vital 
coefficients to be strictly constants. We consider these in turn: 

The coefficient (C). As already stated, we rather expect that 
this coefficient shall diminsh during life in the case of a “typical” 
vertebrate organ such as the liver or pancreas. The factor K of C 
should furthermore be labile with respect to radical changes in type 
of metabolism or to very serious geometrical rearrangements as in 
early foetal life. 

The coefficient (B). This is the vital coefficient related directly 
to the gross metabolism of the cell, and consequently its mode of vari- 
ance will depend on the type of organ. The following general remark, 
however, would seem to hold true for all organs: There should be a 
sudden rise in B coincident with the initiation of actual function, and 
again the initiation of the thyroxine effect should be reflected in an 
increased B. For the rest of life B should be virtually a constant. It 
is quite conceivable that the rapid rise in B should be the best index 
of the inception of function. 

The quantity (a). Two tangible remarks may be about a, and it © 
may well be that they express the operation of the same phenomenon. 
As yet unpublished data suggests that, in the liver, growth hormone 
causes a specific rise in a, and lactogenic hormone likewise but to a 
lesser degree. The mechanical theory of the cell developed by Profes- 
sor Rashevsky and his associates predicts that in the fundamental 
respiratory metabolism of the cell any shift favoring anaerobic gly- 
colysis should increase dN/dt. This effect should be specifically ob- 
servable in a. 

The quantity (Rx/N). This is a coefficient of considerable inter- 
est, for it embodies the effect of the circulatory system and forms the 
basis for a type of experimental analysis which is discussed elsewhere 
(Kreutzer and Morales, in preparation). The circulatory supply 
might be said to be “in gear” when it grows in proportion to the num- 
ber of cells which it feeds and whose wastes it removes. Mathematic- 
ally this amounts to setting R = kN. In this case Rx/N is a true 
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constant as long as the composition of the blood does not change in 
nature (a, constant). There are certain important cases in which 
the circulation runs “out of gear” with the organ. One of these occurs 
in the early stages of development, when the organ cells proliferate 
ahead of blood vessels; a second such circumstance arises following 
the extirpation of part of the organ, leaving the blood supply intact; 
finally, in one classic case at least (Higgins, Mann, and Priestly, 1932) 
it was possible experimentally to enrich the blood supply to an organ 
and observe the result on growth. In all of these cases R > kN with 
time. With the exception of purely temporary disengagements like 
these we believe that R = kN, and T/N is a true constant. 

In closing we should like to reiterate that studies of this sort, 
with the aid of experimental data, may aid in elucidating the general 
nature of population growth. They also are illustrative of methods 
that offer the possibility of drawing inferences regarding individual 
unit behavior from observing the more accessible behavior of the 
population. In the latter capacity, we hope they may serve the solu- 
tion of the “major problem”, namely, to understand the mechanism 
of the single cell. 
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. The standard approximation method used in mathematical biophysics 
is applied to the problem of flow of an incompressible viscous fluid in an 
elastic distensible tube. It is found that the wall of the tube may per- 
form damped transversal harmonic oscillations due to that flow. The 
phenomenon is independent of the viscosity, the latter contributing only 
a damping factor. While, due to rather rough approximations, the 
practical applicability of the equation derived is rather limited, it is sug- 
gested that they may give a clue to the understanding of vibrations of 
the walls of blood vessels which long ago have been suggested as the 
possible source of some hemic murmurs. 


The hydrodynamics of viscous flow still presents a large number 
of difficult theoretical problems, even in relatively simple cases. The 
hydrodynamics of blood flow offers additional difficulties because of 
the non-rigidity of the blood vessels. In any ordinary hydrodynamical 
problem, even for ideal fluids, the shape and size of the tube through 
which the fluid flows enters into the boundary conditions of the sys- 
tem of partial differential equations which describe the flow. In case 
of a distensible tube like a blood vessel, the size and shape of the _ 
boundary are themselves determined by the distribution of velocities 
and pressures which in their turn are determined by the boundary 
conditions. Problems of that nature, not unfamiliar to the mathe- 
matical biophysicist, lead to complicated, usually non-linear, integro- 
differential equations or more general equations in functional deriva- 
tives. Exact solutions of such equations may well be considered as 
hopeless at the present time. It is therefore indicated to use some 
crude approximation method similar to the standard one in mathe- 
matical. biophysics. The results thus obtained represent only very 
rough relations between average quantities, but they give at least 
some insight into the phenomena studied. 

The most important contributions to the hydrodynamics of blood 
tlow have been made in the classical investigations of O. Frank (1899; 
1905; 1920; 1926). Frank’s main interest centers on the problem of 
propagation of pressure waves. The problem of vibrations of the 
elastic walls of the blood vessels is also to some extent treated by 
Frank. Frank goes back to the differential equations of viscous flow, 
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and simplifies them by making some plausible idealized assumptions. 

In this paper we shall attempt to apply to this field the approxi- 
mation method which has been so successfully used in other problems 
of mathematical biophysics. The problem of vibrations of the blood 
vessels can, as we shall see, be treated rather simply by this method, 
leading to expressions for the frequency which is simpler than the 
one obtained by Frank. In subsequent papers we shall treat by a 
similar method different aspects of the problem. The results of this 
paper indicate the possibility of vibrations with frequencies within 
the audible range and thus suggest a possible explanation of some 
hemic murmurs. 

As a preliminary problem we must discuss the relation between 
the radius of the blood vessel and the pressure which its elastic wall 
exerts on a liquid contained inside. 

For simplicity, consider the blood vessel as an elastic tube 
stretched radially along a sufficient length so that its form remains 
essentially cylindrical. 

Let R, be the radius of a blood vessel when it is not stretched 
at all. The actual radius R is usually several times larger. Consider 
a piece of blood vessel of length 1. If under the action of an internal 
pressure p the radius increases by the amount dF, the work done by 
p is equal to 2npRdk . If 6 is the thickness of the wall, then the axial 
cross-section is equal to 6, and let 6 << R,. The relative elongation 
of the circumference when the vessel is stretched from R, to R is 
R—E& 

0 
sume Hooke’s law to hold for rather large relative elongations, we 
have for the elastic force 


= Hence, if E is the elastic modulus of the wall and if we as- 


R—R, 
| ee 


When the radius increases from R to R + dR, the circumference 
increases by 2xdR , and the work done by the force F is equal to 


F=E6 


(1) 


2nE 6d ~ dR. (2) 


0 


Let p- be the elastic pressure due to the force F . When the radius 
Increases by dR, the work done by p, is equal to —2ap.RdR. This 
must be equal to expression (2). Hence 


saree tere Ley =) (3) 
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For R = R, , when the tube is unstretched,» =—0. WhenR=~, 


_E6 

Reo 
In other words, assuming Hooke’s law to hold over any range of 
stretching, equation (4) gives the critical pressure at which the ves- 
sel would burst. Actually, due to deviations from Hooke’s law, the 
vessel may break at much lower pressures. 

For p < p- the equilibrium radius R* of the tube is determined 
from equation (3), in which p is considered as given. 

If a viscous fluid flows through the tube in a stationary state, the 
pressure p drops in the direction of flow. Therefore, the radius R* 
will vary along the axis. The law of this variation will be discussed 
in a subsequent paper. Too rapid a variation of R* in the axial direc- 
tion will result in additional elastic forces due to longitudinal stretch- 
ing of the wall, and the expression for p, will be different from that 
given by equation (3). We shall confine ourselves here only to cases 
when F* varies so little along the axis that at each point equation (3) 
can be applied. Denoting by po; and 2 the hydrodynamical pressures 
in the stationary state at two points, 1 and 2, and by FR," and R,”* the 
corresponding radii, we have 


1 il 1 1 i 6) 
» - ES yi = ==(0; De — E36 aay nas E=Q. 
7 (i =) (ee et 


A periodical deformation along the tube may be represented very 
roughly for very small amplitudes by a set of straight line segments, 


De = De (4) 


---j-- —--j — 
FIGURE 1 


as shown in Figure 1. The quantity 2] may be considered as the wave 
length. Defining R, and R; as indicated in Figure 1, put 
Ae Bee hs Ae= Kz— Ks" , (6) 


and let A, and A, be very small. Then we may neglect possible addi- 
tional elastic forces due to the radially directed parts of the tube wall, 
and still apply equations (3) and (5). For each half wave length of 
the tube we may now write the equations of motion. Denote by p, the 
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density of the wall, and by p, and p» the actual pressures which will 
differ very slightly from the stationary state values po. and Doz, 80 
that 

D=Pu t+ (Ap)i; DP2=Po2 + (A D)e; (7) 
the quantities (A p), and (Ap). being very small. The equations of 
motion now are: 


d? R, 1 1 
Sn Rib et Sales pies ( re! . 


dt? Ro Rk, 
(8) 
d? R. fs 1 1 
ia R, 8p, i ——-=2n lk | ee 
7 ive O pi dt? 2 | ( R, R, ) 
or 
d? R, 1 y | 
6 p1 =p,—E6 ees easy 
dt? Ry R, 
(9) 


CR 1 1 
6 pi ps — Eta 
dt? Fig, «Eee 


Introducing equations (6) and (7) into equations (9), expand- 
ing the right hand sides, preserving only terms linear in 4; and 42, 
and remembering equations (5), we find 


d? A, E6 
6 pi a A, + (A P)s3 
dt? y the: 
if 
dA, E6 410) 
6 pi ee Tc, ee (A D)2. 


Unless the wave length 21 is very large, R,* and R,.* will be prac- 
tically equal. If we consider under those conditions such a deforma- 
tion of the tube which preserves the constancy of the volume, then 

A, = As . (11) 

Introducing this into equations (10), subtracting the second 
equation (10) from the first and putting (1/R,*? + 1/R,"*) = 2/R%, 
where R* is an average value for the whole wave length, we find: 


d? A, 2E 6 
26 pi ie os + (Ap). — (Ap)e. (12) 
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As the tube in region 1 contracts at the rate dA,/dt, the total 
volume of the region varies at a rate 22 R* 1 dA,/dt. This expels a 
total amount of liquid, of density p, 2 2 p R* | dA,/dt from region 1. 
If the total flow through the tube is maintained constant so that at 
the left end (Figure 1) the velocity is constant, this increases the ve- 
locity at the right end of region 1 by — (2 1/R*) (dA,/dt). Hence de- 
noting by v, the average velocity in region 1, and by v the average 

_constant stationary velocity, we have 


et gab 
V34=Vv-— : (13) 
R* dt 
Similarly, we have for region 2 
21dA, Ll dA, iedidh 
V2=v - > OY - — 5; (14) 
Feeaatit R* dt R* dt 
Hence 
dv, lL dA, dve L d? A, 
—=— ,—-—=- = (15) 


dt R* dv dt R* dt 
and therefore, for the derivative of the average velocity v = 
(v, + v2) /2 over the whole wave length we have 


da Ia, 


dt R* de 


We shall now use an approximate form of the Stockes-Navier 
equations for the motion of liquid contained in the region AB (Fig- 
ure 1). 

The.total acceleration of a fluid particle is the sum of the “‘sub- 
stantial” acceleration 6v/édt and of the ‘“‘convective” acceleration 


(16) 


v grad v =4 grad v’. (17) 


The average substantial acceleration of the liquid may be put equal 
to di/dt. The average value of the convection acceleration is 


Vo" ine Wi (Woe as Da 
: Te ea (18) 
l 2l 
Hence the total inertial force acting on the fluid is 


= dv (I? PRS 
akin T+ J. (19) 
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The applied force consists of two components: the force due to the 
pressure gradient, and the force due to friction. The pressure differ- 
ence is p, — 2 and the total force is 


n RY (D1 — Do). (20) 


To compute the force of viscous resistance, we proceed as follows: 
The actual axial velocity varies from the axis to the wall, being 
zero at the latter. The locus of the average velocity 0 is somewhere 
halfway between axis and wall. Hence the average radial gradient of 
velocity is approximately 26/R*. This produces a force per unit sur- 
face equal to 2n0/R*, 7 being the viscosity. The surface along which 


this shearing force acts is 4a R* l. Hence the total force of resistance 


is equal to 
—Suyid. (21) 


Putting the inertial force (19) equal to the sum of expressions (20) 
and (21), we find 


di p(v.2 — v2) 8nld 
Dent em (22) 


Introducing expressions (7) and (16) into equation (22), re- 
membering that v, = v., and that in the stationary state we have 


Sylv 

Re 

where v is the same as in equation (13), so that v —- 3 =v—-—v,= 
(l/R*) (dA,/.dt), we obtain 


Po — Po2 = 5) (23) 


p B d? A; 8 » I? dA, 
ein mal plac pls ee 
R* d# te 
Combining this with equation (12), we obtain 
Pop\ ea 2E6 8nP dA 
(2on+3 nae Bo tg a = af (24) 
Vins dt? ie Fe dt 


This is an equation of damped harmonic oscillations. As may have 
been expected a priori, the role of the viscosity is to add only a damp- 
ing term. With 7 = 0, and with a deformation which preserves the 
constancy of the volume per wave length, the system may formally be 
considered a conservative, since the amount of energy flowing into 
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the tube at one end is equal to the amount flowing out at the other. 
In that case, equation (24) without the damping term can be readily 
obtained by the Lagrangian method. The kinetic energy of the vi- 
brating walls is per wave length 


dA, \? 
T= 203ml ). (25) 
dt 


the bar over R* being omitted because now in a stationary state 
k* = const. The kinetic energy T. of the liquid is, per wave length, 
because of equation (13): 


T2=42R*? 1 p(v,? + V2?) 


Ztv ada, PU-aA, \* 
=7 L p — se ——- ; 
Lie Oe Ne ae 


The potential energy P due to the elastic forces is obtained by remem- 
bering that during the displacement A, the average force acting on 
the displaced part of the wall is 


1 1 
cide AA CHC | cee nate eereer iy fae 
Re eleh asa 


Hence, per wave length: 


1 1 Z 1 
fe ees | iis (aaa oF 
Re | ie a Fe A; KY tA. 


(26) 


27 
2nalH 6 ze 
Pian 
Puiinge 7,— 7, +7. ,and. writing 
Gar. ao oo teh Bi 
bet ae oer ee ee ; (28) 
dtdA, dA dt 


we obtain equation (24) without the damping term. 


A similar result may be obtained for 7 = 0 by using exact hydro- 
dynamical equations for infinitesimal sinusoidal deformation of the 
tube. This will be shown in a subsequent paper. 


Neglecting the damping term, equation (24) gives a frequency 
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2B 6 
r=] ————.. (29) 
23 p RY +E pR 


With R* <1 cm, 6 © 0.1 cm; p © pp. © 1 gmecm® and/~™ 1 cm, this 


simplifies to 
2E 6 
= —_—— 30 
v4 f 2 R R* ? ( ) 


and with E © 10°, it gives vy © 102, that is, within audible range. 

It must be remarked that for 7 = 0 there may be two values of 
R* for the stationary state when the tube remains cylindrical. In 
that case, p is constant along the tube, and, if p, denotes the constant 
pressure outside of the input, we have 


pv 

me 
If the total flow Q of liquid through a cross-section of the tube is con- 
stant, then v = Q/z p R* and therefore 


DP = Do — (31) 


Q? 
| ad |) aperomemmentnr sel (32) 
2 7 p R* 
Equations (5) now become, since Po: = Po2 = p: 
E 6 Q? E6 
iL i charrage tne pumarme sy 7 =0. (33) 


For reasons discussed in connection with equation (4), p — Hé/R, 
<0. Therefore, as can be seen graphically, equation (33) has two 
positive roots, R* and R*” > R*’. Since the left side of equation (33) 
represents the total force acting radially upon the tube, it is also seen 
that R*” is stable, while R* is unstable for such deformations of the 
tube, which leave it a cylinder. For wave-like deformations, consid- 
ered above, both may be stable. Physically, however, oscillations are 
not likely to occur around R”. 

We do not wish to suggest that the hemic murmurs are neces- 
sarily due to such transversal vibrations of the walls of blood vessels. 
The experimental studies of T. Weber (1855) suggest the possibility 
of their being due to longitudinal elastic vibrations of the blood: ves- 
sels. According to Weber, the mechanism may be somewhat similar 
_ to the action of a uniformly moving bow upon a violin string. Possi- 

bilities of relaxation oscillations must be definitely considered in this 
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connection, and further theoretical studies are indicated, of which 
this is only the first step. 

This work has been aided in part by a grant from the Dr. Wal- 
lace C. and Clara A. Abbott Memorial Fund of the University of 
Chicago. 
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In an elastic distensible tube, like a blood vessel, the radius is de- 
termined by the equality of the hydrostatic pressure and the elastic 
forces. If a viscous fluid flows through such a tube, there is a pressure 
drop along the line of flow. This results in a variation of the radius 
of the tube along the axis. An approximate expression, valid within a 
limited range of values, is derived for the radius of the tube as a func- 
tion of the distance along the axis. Another approximate expression is 
derived for the relation between pressure drop and total flow in such a 
case. For sufficiently high rates of flow the pressure drop does not vary 
linearly with the flow, as in the usual Poiseuille’s law, but more rapidly. 


In a previous paper (Rashevsky, 1945) we investigated the vi- 
brations of the walls of an elastic tube through which a viscous fluid 
is flowing. We found that the frequency of those oscillations depended 
on the radius of the tube. But in case of a viscous liquid, that radius 
itself varies along the axis of the tube because it is partly determined 
by the hydrostatic pressure and the latter drops along the line of flow. 
In this paper we shall investigate the law of variation of the radius 
of the tube for the case of a stationary flow. The exact solution of 
the problem presents very great difficulties. It amounts to solving the 
Navier-Stokes equations for an arbitrarily shaped axially symmetric 
tube, and then determining the arbitrary function, which describes 
the shape, from considerations of mechanical equilibrium. We shall 
therefore suggest here a very crude approximation which will give 
us a general idea about the variation of the radius, at least within a 
limited range of values. 

Introducing cylindrical coordinates, r, z, and ¢, we first of all 
restrict ourselves to axially symmetric flows, so that ¢ and its deriv- 
atives drop out from ali equations. Further, we shall restrict our- 
selves to the case that the radius R of the tube varies very slowly 
with z. 

Our approximation amounts essentially to two things: first we 
shall neglect the very small radial component v, of the velocity v 
everywhere, and second we shall assume that at each cross-section of 
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the tube the radial variation of the axial component v is parabolic, 
like in an ordinary Poiseuille flow. We shall omit the subscript z and 
write v for v,. For a stationary state, when all time derivatives are 
zero, the Navier-Stokes system of equations reduces now to a single 


equation: 
Op ov ov. Lda ee 
-—= 74 ++ <-— + ). (1) 


where p denotes the pressure, p the density of the fluid, and 7 the 
viscosity. 

Denoting by v, the velocity at a point of the axis, and remember- 
ing that everywhere for r = Rk, v = 0, we have 


oe 
=— Vo = a e 2 
v= = (2) 


The total flow Q, which in a stationary state is the same for all 
cross-sections, is given by 


R 
Q=22p [ eran. (3) 
0 
Introducing expression (2) into expression (3), we find: 
2Q 
ap R? 


As a first approximation, we shall consider that R varies with 
z so slowly that the resulting variations of v with respect to z may 
be neglected, in other words, that we can put dv/dz = 0?v/dz? = 0. 
Introducing now expressions (2) and (4) into equation (1), we find: 


Q=42pv R?, or y= 


(4) 


dp 87 
dz ap Rt (5) 
which gives 
87Q [7% dz 
D= Fa : 
mp . R*(z) (6) 


Denoting by 6 the thickness of the tube wall, by Z its modulus of 
elasticity, and by R, the radius of the tube when completely un- 
stretched, we have as the condition of mechanical equilibrium at 
every point [loc. cit. equation (5)]: 
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a i: 
o—E6{. ——-— | =0. (7) 
he Ee 
Introducing expression (6) into (7), we have 
87Q er A E56 E 6 
Do — — + = (0. (8) 
np F R+ (z) Ro R(z) : 
Differentiating this with respect to z, rearranging and putting 
Eéap 
=A, (9) 
87Q 
we obtain 
dR 1 
A—~=- —. (10) 
dz R? 
Integrated, this gives, denoting by B a constant of integration, 
R Zz 
—=B--, (11) 
3 A 
or 
= 32 
R= 5 8 Pee cpt (12) 


Because of the original assumption which we made about the slow 
variation of R with z, this expression holds only within such ranges 
of z for which z/A <<B. In this case, we may expand the right side 
of equation (12) preserving only linear terms in z. Putting 


1 87Q 
——— ———_ EC, (18) 
AB EéapB 
we obtain 
a 
R= (ap)"( 1-22 . (14) 


The quantity (3B)? is the radius for z= 0, or forp—=p,. It 
is determined as the root of equation (7), into which we put p = py. 
Denoting that root by R;, we have 

E 6, Rk. 247Q 
R; = ————__; B=—-; C=———. (15) 
E 6 — Do Ro 3 E bap 
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Substituting expression (14) into equation (6), and making use 
of the second relation (15), we find after integration: 


87Q 1 1 
p-p=—— [——-¢ )- (16) 
np hy C(t 


Expanding (1 — C z/8)*, keeping only the linear terms, and using 
expression (15) we finally find 


87 Qz Pe 
Do Doin ao 5 
eee eae 
Eéapk3 


For very small values of Q the relation between p, — p and either Q 
or z is linear, as in Poiseuille’s law. For larger values, however, it 
is not so. If R; were constant, then p. — p would vary more rapidly 
than Q or z. But R; is itself a function of p,. We still may keep 
»,, and therefore also RF; , constant, and then see that the total flow 
varies more slowly with the total drop of pressure than in Poiseuille’s 
law. If we, however, keep p constant and vary p, the situation be- 
comes much more complicated. The flow then is not a function of 
the pressure difference only, but also of the absolute value of the 
pressure. The actual relation is obtained by substituting expression 
(15) for R; into equation (17). 

We may obtain a better approximation by dropping the restric- 
tion dv/dz = 0, but still retaining 0?v/déz? = 0. We now introduce 
into equation (1) an average value dv/dz, defined in the following 
way: 


Let 
dR 

z)=—. 

f(z) Pm (18) 
Then we have 

ov soav . 
on ae oa Cd (19) 
oz OR 


Introducing expression (4) into expression (2) and the latter 
into (19) and averaging over the whole radius R, we find 


w_8 | 
2 Qf (2) ['re- ue fe= 4QF (2) on: SR 


np R® 7 p te* 3apR> 


Introducing this into equation (1) and remembering that the 
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average value 0 of v is, from equations (2) and (4), equal to 
0 = 4Q/a p R?, we now find 
dp = 16Q°f(z) 87Q 


dz 872pR> apR 


Integrated, this gives 


16 @ (*7(z)dz 8,Q [* dz 
Pp=P. + = é 22 
32° p} R(z) 7p , 2 (z) ae 
Introducing this into equation (7), differentiating with respect 


to z, and rearranging, we obtain 


= 1 SD es 1 


(21) 


3722p R>  R? 


Integrating, and denoting by C’ a constant of integration, we 
have 


=, (23) 


1 16 Q? 87Q 
—foh = lor k= GC’ — 
3 3 7 p Lp 


Bhs (24) 


If we omit the second term of the left side, we obtain again the 
expression (11). If that term is, however, retained, we find that R 
is a two-valued function of z. The left side is positively infinite for 
Rk = 0 and for R = ~, and has a minimum. Below that minimum 
value of R, no real values of z exist. In the neighborhood of the 
minimum, however, the solution is not valid since here we cannot 
neglect 0?v/dz?. The existence of two values of F& is of the same 
nature as a similar phenomenon which happens for 7 — 0 and which 
has been discussed in the previous paper (Rashevsky, 1945). Only 
the larger value of R is stable for 7 — 0, and this is likely to remain 
true for 7 > 0. Sufficiently far away from the minimum, and keep- 
ing on the larger branch of R’s, we again obtain equation (11) with 
all its consequences, for the term in R* prevails over the log R term 
in equation (24). 

lf we wish to obtain a still better approximation, we may drop 
the assumption 0?v/dz2 = 0, and introduce into equation (1) an 
average value 0?v/dz? obtained in a similar way as 0v/0z. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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In continuation of previous studies, inequalities between different 
parameters of the brain are derived which determine whether an indi- 
vidual prefers in general visual patterns consisting of a relatively small 
number of relatively strongly excited elements, or such patterns which 
consist of a very large number of relatively weakly excited elements. 
As has been discussed in a previous publication, the first type of pattern 
is usually represented by artificial human-made designs, whereas the 
second type of pattern is formed predominantly in natural landscapes 
and sceneries. Thus the inequalities established in this paper give us 
the biophysical conditions which determine an individual’s preference for 
either artificial designs or for landscapes and other natural objects. 


In preceding publications (Rashevsky, 1938b, 1940, 1942a, b; 
Rashevsky and Brown, 19-44a, b), we have developed a biophysical 
theory of visual aesthetics based on the assumption that the indi- 
vidually felt aesthetic value of a visual pattern is proportional to 
the total central excitation produced by a given pattern. Applied to 
relatively simple artificial designs, the theory leads to predictions 
which are confirmed by experiment (Rashevsky and Brown, 1944a, 
b). It has been pointed out (Rashevsky, 1942a) that two distinct 
cases must be considered: one is the case of relatively few elements 
excited sufficiently strongly so that the thresholds of the inhibitory 
cross-connecting fibers may be neglected; the other is the case of a 
very large number of elements excited sufficiently weakly so that the 
inhibitory cross-connecting fibers remain unexcited. Artificial man- 
made designs are typical representatives of the first case, while land- 
scapes are characteristic of the second (Rashevsky, 1942a). How- 
ever, mixed cases do exist, and their occurrence is not infrequent. 
Thus in a landscape the aesthetic effect of the very large number of 
small differences in shades and hues may be just as important as the 
general design, or what the artist refers to as “composition”. It is 
quite legitimate and theoretically very interesting to inquire whether 
a difference in the biophysical structure or biophysical parameters 
of the central nervous system may result in a preference by an in- 
dividual towards one or another of the two types of patterns men- 
tioned above. The problem has more than a purely theoretical inter- 
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est since actually different individuals display such preferences. One 
artist may pay much more attention to composition, another to 
color effects. An individual may definitely prefer a tree to a 
flower or vice versa. Now in a tree the geometric “design” plays a 
relatively subordinate role, while the effect of the very large number 
of different shades of green or other colors is more outstanding. Even 
the general outline of a tree consists of a very irregular line, with very 
many not too sharply defined elements. In a flower, on the contrary, 
the element of a geometric design is very much more pronounced. 

Let-us investigate the problem from the point of view of the 
previously developed theory. We shall use notations employed previ- 
ously (Rashevsky, 1938a, Chapter xxii, referred to as loc. cit.; 1942a) 
and shall presuppose familiarity with the previous publications. 

It is of course impossible to attack at present the problem in its 
generality. However, even some approximate, abstract, theoretical 
considerations are of interest. We shall therefore treat the case when 
all elements of a pattern are stimulated with the same intensity and 
all neural elements involved have the same biophysical constants. 

For a “natural” pattern consisting of a large number of weakly 
excited elements we must always have (Rashevsky, 1938a; 1942a) 


S; # hy 


. (1) 
bos AEF 


If m denotes the total number of peripheral nerve fibers involved or, 
with our assumptions, the total number of elements, then the total 
central excitation is given by multiplying by n the first term of the 
right side of equation (53) of loc. cit., since, because of expression 
(1), H;, and hence the second term, are zero. Doing this and intro- 
ducing expressions (10) and (61) of loc. cit., we find for the total 
central excitation 

E=nPa,lI,(S; — hi). (2) 


This expression increases with x. The consequences of this have 
been discussed previously (Rashevsky, 1942a). For a fixed it in- 
creases with S,. But the latter is subject to restriction (1). Hence 
for a fixed n the largest possible value of E is given by 


5, ak: (3) 
Consider now the case when 
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SO that the inhibitory cross fibers are excited. Then the total central 
excitation EH’ is obtained by multiplying by n the right side of ex- 
pression (60) of loc. cit. This gives 


E’ =n{P a, 1,(S: —'kh) 


—(nh— 1) Qa, 1, Pa, 1, (S, — h,) — hy]} . Ae? 
For a fixed S,, the quantity H’ has a maximum for 
rim ey Por a) ae de LP oo, J, (Sy — hy) ha] 
a 2 Q as Iy[P oy Iy(S, — ha) — he] Daatty 
and that maximum value is equal to 
ae os, —h,) 7 Oe, af, (8;— kh.) — hy]}* es 


4Qa;1;[P 0 1, (8; =aitg) hs] 


The value of E’,, increases with S,, and at first glance may be made 
infinite by taking a sufficiently large S,, Actually this is not so. Ex- 
pression (7) is derived on the assumption of a linear relation between 
E, and S,, and such a relation holds only for not too large values of 
S,. As S, increases, #, tends to its maximum values I,/6,, 0, being 
the refractory period. For very large values we therefore must sub- 
stitute I,/6,, for a, I,(S, — h,) into equation (7), and thus obtain 
{PI,/6, + Qa; 1,[P 14/6, — he] }? 
a (8) 
At nile (i 1 bee | 


Equation (8) holds still subject to the assumption that while the 
excitatory fibers have reached their saturation values of excitation, 
the inhibitory fibers are still within the linear range. If the latter 
are also near their saturation, we have, instead of equation (8), 


[PL /t + Q1:/6] 
Set 14/85 >) 


Finally the case must be considered when the inhibitory fibers 
reach their saturation, while the excitatory are still in the linear 
range. We then have 


[P o, 1,(81 — K) + Q1;/63]? 
2 
4Q 1/95 


In this case again E’,, increases with S,, and its maximum value is 


(9) 


(10) 
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reached when E, becomes equal to J,/6,. Then equation (10) reduces 
to equation (9). 

If E’, is approximately equal to E,,, then the individual has no 
particular preference for either “natural” or “artificial” patterns. 
If E'm >> Em, the individual prefers artificial designs; if E’m << 
E,,, the person prefers “natural” patterns. Since the maximum pos- 
sible value n,, of n in equation (3) is a constant for a given indi- 
vidual, the above inequalities give us the relations we are looking for. 
In case of equation (8) the condition for preference of geometric de- 
signs is given by 


{P 1,/6, + Q a3 Is[P 11/81 — hs] }? 


o> Naks. (11) 
4Q a; I,(P 1/61 =o hs] 
In case of equation (9) we have 
[P 1/0, + Q I;/63]? 
> hr hts. (12) 
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In both cases low threshold h; of inhibitory fibers favors prefer- 
ence for geometric designs. For a fixed h; a large nm», which may be 
interpreted as a more developed central nervous system, favors pref- 
erence for “natural” patterns. For very large J; inequality (12) be- 
comes 


Ts 


4 63 


>> Mm hs, (13) 


indicating that very strong inhibitory fibers favor preference for 
geometric designs. 

For a fixed value of J,/6,, the left side of inequality (12) has a 
minimum for 


| Rene? aad & 
6,5 QA 
the value of that maximum being P J,/6, . 
If 
PI,/6, > Ton hs, (15) ; 


then for values of I;/@; which are sufficiently smaller or larger than 
those given by equation (14), inequality (12) will be always satis- 
fied, the individual prefering geometric designs. Even if inequality 
(15) does not hold, inequality (12) holds for sufficiently small or suf- 
ficiently large values of I;/6,, since the left side of equation (12) is 
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infinite for I;/6; = 0 and for I;/6; = ©. But in general, P/Q ~ 1. 
Hence equation (14) implies J;/@; © I,/@,. In other words, an in- 
dividual with either too weak or too strong inhibitory neurons will 
prefer geometric designs. For the preference of “natural” patterns 
or for lack of definite preference, the saturation excitations of ex- 
citatory and inhibitory fibers must be approximately equal. 
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